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Abstract 


We tested a new model of CMOS detector manufactured by the Gpixel Inc, for potential space astronomical 
application. In laboratory, we obtain some bias images under the typical application environment. In these bias 
images, clear random row noise pattern is observed. The row noise also contains some characteristic spatial 
frequencies. We quantitatively estimated the impact of this feature to photometric measurements, by making 
simulated images. We compared different bias noise types under strict parameter control. The result shows the row 
noise will significantly deteriorate the photometric accuracy. It effectively increases the readout noise by a factor of 
2-10. However, if it is properly removed, the image quality and photometric accuracy will be significantly 


improved. 
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1. Introduction 


Since the time of digitalization, CCD (charge coupled 
device) detectors have played an essential role in astronomy 
(Mackay 1986; Janesick 2001; Howell 2006; Lesser 2015). But 
after the invention of CMOS (complementary metal oxide 
semiconductor) image detector (Fossum 1997), CMOS has 
attracted more and more attentions in astronomy due to its 
advantages (Qiu et al. 2013; Wang et al. 2014; Mendikoa et al. 
2016; Jorden et al. 2017; Shugarov 2020; Ardilanov et al. 
2021; Greffe et al. 2022; Liu et al. 2022), especially in time 
domain area (Pratlong et al. 2016; Niu et al. 2022; Song et al. 
2022). It also shows the potential to replace CCD (see e.g., 
Gach et al. 2022). 

During the development of the Chinese Space Station 
Telescope (CSST), both CCD and CMOS were considered as 
possible options for the CSST Survey Camera (CSC). The 
CSC is designed for large area optical survey. It uses 30 
detectors to form a giant 2.3 billion pixel array to cover a 
field of view of more than one square degree. The optical 
survey will cover wavelengths ranging from 250 to 1000 nm. 
The detectors are required to have high quantum efficiency, 
low readout noise, low dark current and large dynamic range 
to fit for the proposed wide field blind survey (Zhan 2011, 
2018, 2021). In order to test whether the state-of-art CMOS 
detectors can satisfy the scientific requirement, we perform 
series of tests in the laboratory at the National Astronomical 
Observatories, Chinese Academy of Sciences, in Beijing. The 
detector we tested is an HR9O9OBSI scientific CMOS 


(sCMOS) detector? manufactured by the Gpixel Inc.° This 
model is under active development in last few years. The 
Gpixel Inc. also produces some other science level CMOS 
detectors for (potential) astronomical applications (Shugarov 
2020; Gill et al. 2022). 

The readout noise is one of the most important parameters to 
evaluate the performance of a detector. The bias noise properties 
do have subtle consequences to final photometry or spectroscopy 
quality (Jiménez Mufioz et al. 2021). Only a few studies were done 
on CMOS detectors in astronomy field (Basden 2015; Basden & 
Morris 2016; Greffe et al. 2022), probably due to limited 
application cases. More can be found in non-astronomy scientific 
applications (Diekmann et al. 2017; Matsiaka & Plakhotnik 2021). 
In this paper, we will study the bias noise properties of this detector 
and its photometric performance with simulated images. Since this 
model is still under development, the content of this paper only 
reflects the status before the end of 2022. 


2. The Noise Properties of CMOS Bias Image 
2.1. The CMOS Detector and the Bias Images 


The bias images were taken from an HR9O90BSI sCMOS 
detector (serial number HR9O90BSI-AIOLV-SV-22001JC- 


> The detector we tested is a customized one. The standard version of this 


model is officially named as “GSENSE1081BSI.” A Chinese exoplanet 
searching project called “Earth 2.0” (Song et al. 2022) uses the standard 
version. 


© official website: https: //www.gpixel.com/. 
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Figure 1. The setup of the detector test bench. 


0027). This one is designed to optimize the performance in 
visible bands. The pixel physical size is 10 wm x 10 um. The 
detector supports on-chip 16 bit column-parallel analog-to- 
digital converter (ADC), with up to four low voltage 
differential signalling (LVDS) channels for image output. 
The detector uses rolling shutter. The supported minimal 
exposure time is about 0.6 ms, which is used for taking bias 
images. The readout time is about 3 s when using two LVDS 
channels at frequency of 250MHz. The raw image has 
8976 x 9222 pixels. On both left and right sides, there are 34 
columns at the edge of array designed as “electrical black” 
(EB) region. The readings of these pixels only come from 
electronics. In principle, they can be used to calculate bias level 
as well (similar to overscan region in CCD case). However in 
practice we found the values in EB region were exceptionally 
large, and therefore not useful for our analysis. In this paper, 
we only use the central 8900 x 9120 pixels (“active” region). 

We prepared a test bench (Figure 1) for the test of CSC 
detectors in an ISO 7 clean room. It can satisfy multiple test 
requirements (e.g., pixel response at different wavelengths). 
The optical system was kept in a light tight box, which was 
installed on a vibration controlled optical table. All devices 
emitting light (e.g., running indicator) were carefully separated. 
Light intensity monitor system was also used when necessary. 


It makes sure the bias and dark images are not contaminated by 
any external light source. The detector was installed in a 
vacuum dewar, cooled with liquid nitrogen. The temperature of 
the detector was kept at about 188 K, with typical temperature 
fluctuation no more than 1 K, similar to the designed on-board 
environment. 

The bias pattern is not very sensitive to the temperature. The 
variation of temperature has negligible contribution to the bias 
noise. The backend test board was also provided by the Gpixel 
Inc. (serial number 008). The test was performed in full dark 
environment to avoid external light contamination to the bias 
images. Given the typical dark current of 0.02 e~ /s/pixel and 
high readout frequency, the dark current contribution during 
readout is negligible. In this paper, we use 30 bias images that 
were continuously taken on 2022 September 5. We note that 
other data sets show similar properties. 

In this paper we will use ADU for simplicity. The gain value 
(in astronomical definition) is about 1.3e /ADU (measured 
from photon transfer curve’). Readers can convert all the 
results with this gain if they prefer unit of electrons. 


7 Tn the CMOS detector, each pixel has its own gain value. The traditional 
photon transfer curve technique only gives a typical value. 
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Figure 2. Left: master bias. Middle: individual bias image. Right: noise residual image (middle-left). 


2.2. Noise Residual Images and Row Noise 


We first stacked the 30 bias images to make a master bias. 
The stacking was performed pixel by pixel. In each pixel, we 
sorted the 30 values, excluded the top two and bottom two 
values to avoid any possible cosmic ray contamination, and 
then took the average of the rest. We note if we do not exclude 
any data or exclude more data, the result does not significantly 
change. On the master bias, we identified a pixel as “bad” if its 
value exceeded the 70 threshold of typical (median) bias value. 
There are totally 52565 pixels marked as “bad,” including two 
rows. Most defected pixels are excluded in this way. The 
master bias has a strong feature of horizontal and vertical lines 
as well as square blocks. A similar feature was also observed in 
the other CMOS model with large arrays (e.g., Wang et al. 
2014; Karpov et al. 2020). 

The master bias was then subtracted from each individual 
bias image to remove the fixed pattern, leaving only the noise 
fluctuation. Hereafter we call such a fixed-pattern-corrected 
bias image as “noise residual image.” On the right panel of 
Figure 2 there shows one example. The horizontal stripes can 
be clearly observed. We calculated the average value and the 
standard deviation in each row or column on the noise 
residual image. The bad pixels were excluded from the 
calculation. We merge all the row and column values from 
the 30 noise residual images together to make statistical plots 
(Figures 3 and 4). 

The result shows that the row-to-row variation 
(o = 3.59 ADU) is much larger than column-to-column varia- 
tion (ø = 0.05 ADU). It means the bias images contain spatially 
correlated noise, which is visually represented as horizontal 
stripes in noise residual image. The standard deviation in 
each column is comparable with the pixel value noise 
(which is based on Gaussian fitting of the pixel histogram, 
o = 4.46 ADU) and nominal “readout noise” (which is from the 
standard deviation of noise residual image, o = 4.77 ADU). 
Therefore the row-to-row fluctuation is the dominant 


component of the total readout noise. The standard deviation 
in each row is actually quite low (typically ø = 2.84 ADU), 
implying that the detector may have the potential to decrease 
the overall readout noise to less than 3 ADU. 

We notice that the noise distribution is always non-Gaussian. 
The pixel-wise noise distribution is known to be somehow non- 
Gaussian. An internal technical report from the Gpixel Inc. 
(private communication) shows that the histogram of single 
pixel noise values has a long tail at the high end. A similar 
behavior was also reported in other CMOS models (e.g., Ishida 
et al. 2018). 

We calculated the pixel-to-pixel cross-correlation coefficient 
as a function of spatial shift (both x and y direction). The result 
is displayed in Figure 5. The result clearly shows the 
correlation in each row (i.e., x direction) is strong (>0.5). 
There are also many horizontal lines with weak but statistically 
significant correlations, positively or negatively. Hereafter in 
this paper we call such a feature as row noise (RN). We note 
that the pattern also shows some periodicity in the y direction, 
implying that the noise also contains characteristic frequency 
components. 

In order to confirm the existence of the characteristic spatial 
frequencies in the y direction, we calculated the spatial power 
spectra of the bias noise in Figure 6. The power spectrum is 
defined as the squared absolute amplitude of the Fourier 
transform component at given spatial frequency. We calculated 
the power spectrum in each row (or column) and derived the 
stacked spectrum with 20% trimmed mean value for all rows 
(or columns) in each individual bias image. Then we obtained 
30 spectra (illustrated in gray lines) and derived the average of 
them (red line). In this way, we can reject any potential outliers 
(due to bad pixels) to get a robust estimation of the intrinsic 
power spectrum. 

We find that the noise in each row (top panel of Figure 6) has 
a significant low frequency component, consistent with “1/f 
noise” (see e.g., Tian & El Gamal 2000). The noise in each 
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Figure 3. Top right: pixel value histogram of all 30 noise residual images (fixed-pattern-corrected bias images). Bottom left: average value versus standard deviation 
in each row (the Pearson correlation coefficient is displayed). Top left: histogram of average (in each row) values. Bottom right: histogram of standard deviation (in 


each row) values. The red dashed line is Gaussian fitting. 


column (bottom panel of Figure 6) has flat spectrum except that 
there are five characteristic frequencies (0.048, 0.133, 0,181, 
0.229, 0.410, in units of pixel”'). It may be caused by the 
power supply fluctuation, which is quite common in sCMOS 
detectors using column parallel readout architecture (i.e., each 
row is read in sequence, see Mikkonen 2016; Gilroy 2020). 
The readout could be affected by internal or external electrical 
circuits which has characteristic frequencies. However, we note 
that the white noise component is still the dominant noise 
source (the observed white noise component is equivalent to a 
noise source with rms of <4 ADU). 


2.3. Image Correction to Remove Row Noise 


The row noise, or more generally stripe noise, is commonly 
observed in images in remote sensing and/or taken by CMOS 
detectors (e.g., Simpson et al. 1995; Rakwatin et al. 2007; 
Gilroy 2020; Wei et al. 2020). There are many developed 
algorithms to remove the stripe noise, for example, histogram 
matching with facet filter (Rakwatin et al. 2007), midway 
histogram equalization (Tendero et al. 2010), low-rank-based 
single-image decomposition (Chang et al. 2016), group sparsity 
based regularization model (Chen et al. 2017), wavelet deep 
neural network (Guan et al. 2019), etc. It is also observed in 
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Figure 4. Same as Figure 3 but in columns. 


some astronomical data (e.g., Schlawin et al. 2020, for JWST 
NIRCam in time-series). The problem is not new and there are 
many possible solutions. Therefore it will be useful to investigate 
the properties of row noise corrected images as well. 

The origin of row noise in CMOS is complicated 
(Gilroy 2020). In this paper, only additive (or external signal 
independent) row noise is discussed (as only bias images are 
involved). Our focus in this paper is not the correction method 
itself either. Therefore we use a very simple and 


$ Actually, in another imaging test, we find there is row direction crosstalk, 


i.e., signal dependent row noise, when there are many saturated pixels in a row. 
But it only appears when there is significantly saturated, i.e., very bright, target. 
It does not affect the main conclusion of this paper, where we care more about 
faint objects with low signal-to-noise ratios. If we add the crosstalk in our 
simulation (presented in next section), the result will never be better. 


straightforward method: subtracting the averaged empty pixel 
value in each row. 

To start, it is important to define the empty pixels, i.e., the 
image area not exposed to external light’. In case of 
HR9O90BSI, the EB region (34 columns on each side) could 
be used for this purpose, after proper development in future 
(suggested by the Gpixel Inc.). Or alternatively, it is possible to 
mask some active pixels with opaque material to make an 


Image area exposed to homogeneously distributed background light but not 
affected by astronomical objects may also be used for this purpose. However, 
in real application, such exercise has limitations: e.g., it is difficult for large 
extended objects (e.g., Magellanic Clouds) or over-crowded region (e.g., star 
cluster) where hardly any pixel is not contaminated by targets. In this paper, we 
do not discuss the algorithms to overcome such problems. Nevertheless, if 
used, it should be tested in real applications. On the other hand, our solution is 
the simplest way to avoid such problem. 
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optically blind area at the edge of image. This option was 
proposed by the Gpixel Inc. (private communication) as well. 
In this paper, we assume the former option will be available. So 
we test the feasibility of this algorithm here. We use the 
leftmost 34 and rightmost 34 pixels of the active region (to 
mimic the size of EB) in each row to calculate the row average. 
We subtracted the average value of the 68 pixels in each row of 
the noise residual image. A similar algorithm is often used in 
CCD image processing, in which the reference empty area is 
from overscan instead of physical pixels (e.g., Wolf et al. 
2018). In Figure 7 the noise residual image before and after the 
correction are compared. Hereafter, we call these corrected 
images as “RN-corr” images. 

The row noise correction significantly reduces the readout 
noise (from 4.77 ADU to 2.89 ADU). The horizontal stripes 
become much less prominent visually. Similar power spectrum 
and correlation coefficient analysis were performed on RN-corr 
images for comparison (see Figures 8 and 9). 

Although there is still some small pixel-to-pixel correlation 
residual left in each row, the amplitude is greatly suppressed 
(Figure 8). The strength of the power spectrum in vertical 
direction (~6 x 104, Figure 9) is about one third of the original 
one (~2 x 10°, Figure 6), becoming comparable with the 
horizontal one. Though as side effect, the power spectrum is no 
longer flat. Some characteristic frequencies are still there (e.g., 
at 0.048 pixel~'), but their amplitudes are also greatly reduced. 
All these results demonstrate that the correction is very helpful. 
In next section, we will include the corrected images into our 
analysis as well. 


3. The Impact of Row Noise to Photometry 
3.1. The Simulated Noise Residual Images 


Here we use simulated images with galaxies and stars to 
better quantify the effect of row noise in astronomical 
application. We always simulate master bias corrected images 
(i.e., object added to the noise residual image). 

For RN images, we use cutouts from original CMOS 
noise residual images. We randomly pick sub-regions of 
196 x 196 pixels. Each random block must contain no more 
than 10 bad pixels and all bad pixels were set to zero. Random 
blocks should never overlap with each other so that the same 
noise pattern would not appear twice. There are totally 30446 
random cutouts taken from 30 images. We also made cutouts 
for RN-corr images at the same positions for controlled 
comparison. 

For comparison, we also simulated pure Gaussian noise 
images. We made two sets of noise residual images with 
o =4.77 ADU (Gauss-high) and o = 2.89 ADU (Gauss-low) 
respectively. They were used to compare with RN and RN-corr 
images respectively, in sense of nominal readout noise. We do 
not only observe the effect of overall readout noise reduction, 
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Figure 5. Pixel value correlation coefficient as a function of spatial shift. 


but also the effect of RN itself. Figure 10 compares the 
different noise residual types. 


3.2. Real Galaxy Models 


Real galaxies were added to the simulated noise residual 
images. The image cutouts of real galaxies were directly taken 
from the website of the Legacy Survey'°(Dey et al. 2019). 15 
galaxies were randomly picked. For simplicity and better 
signal-to-noise ratio, all of the three bands were merged 
together. Source detection and segmentation were performed to 
remove all nearby neighboring objects (except one galaxy 
a = 187.5740, 6 = 73.0354 where a possible foreground star 
was kept to mimic the effect of possible confusion in real 
pipeline photometry). All pixels out of the target galaxy 
segment were set to zero. The cutout images were then 
normalized to unity total flux. Figure 11 shows all the galaxies 
we used. 

When adding the galaxy model to noise residual cutout, we 
randomly picked one model from the 15 models, multiplied by 
arandom flux scale between 0.1 and 10. We use the same set of 
models for each type of noise residual image, to make sure that 
the difference only comes from noise residual. The left column 
of Figure 10 shows one example of such an image set. 


3.3. Sérsic Models 


Sérsic model is a good way to parameterize the light 
distribution of elliptical galaxies, often used in simplified 
simulations (e.g., Holzschuh et al. 2022). We use python 


i https: / /www.legacysurvey.org/ 
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Figure 6. Spatial power spectra of bias noise in rows (top) and columns (bottom). Grey lines are from individual bias images. Red lines are the average of them. In 


bottom panel, the values of five characteristic frequencies are noted. 


package astropy to make 2D Sérsic model images, with random 
parameters listed in Table 1. Some examples are shown in 
Figure 12. Since the Sérsic model extends to infinite radius, we 
made a segment (consistent with an elliptical aperture) which 
included 90% of the total flux for further photometric 
measurement. Different from real galaxy models, Sérsic models 
have larger variability in galaxy size, probably more repre- 
sentative for field galaxies in a blind survey. The middle 
column of Figure 10 shows one example for different noise 
residual background. 


3.4. Star Models 


To simulate stars, we use 2D Gaussian model with 
o = 1.13 pixels (FWHM 2.66 pixels), which is the designed 
optical resolution (dynamic) of CSC. The model values were 
first assigned to a 0.05 pixel resolution sub-grid and then 
summed together to recover the true pixel value. Random 
amplitudes (between 10 and 10°?) and x/y center offsets 
(between —0.5 and 0.5) were used for each simulated star. 
Some examples were displayed in Figure 13. A circular 
aperture with radius of 4 pixels was used as segment for further 
photometric measurement. The right column of Figure 10 
shows one example for different noise residual background. 


3.5. The Photometry and Basic Parameters 


We used the python package photutils v1.5 to measure 
several basic photometric parameters: position, elongation, 
position angle of the major axis, segment flux, Kron flux and 
circular aperture flux. 

The position measurement is very important in astrometry, 
i.e., determining the actual position of the object. Many other 
more complicated photometric parameters are sensitive to 
position measurement accuracy. The elongation is defined as 
the ratio between the length of major axis and the length of 
minor axis. It is often used as an indicator of the inclination 
angle of spiral galaxy. The position angle of the major axis is 
defined as the angle between the major axis and y-axis 
direction. It is also an important parameter to describe how the 
galaxy is positioned. Both the elongation and position angle are 
important morphological parameters, especially useful in weak 
lensing shear measurement (although the actual algorithm is 
more complicated). 

The segment flux is the sum of pixel values within the 
segment. It is a model-independent flux measurement, which is 
believed to be a robust estimation of the total flux. The Kron 
flux (Kron 1980) is another parameter to calculate the total flux. 
It sums up the flux within an elliptical aperture (i.e., assuming 
the galaxy shape can be represented by an elliptical). The 
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parameters of the ellipse are determined from the first-order and 
second-order moment values. We use the default parameter in 
photutils. Comparing with segment flux, Kron flux is more 
sensitive to position and morphological measurements, but 
more reasonable in sense of measuring the total flux of a 


(generally symmetric) galaxy. It is widely used in many optical 8 

surveys (e.g., Guo et al. 2013; Chambers et al. 2016). The 2 08 
circular aperture flux was also calculated, which was preferred Š 0.6 
for stars. We always used the circular aperture with radius of 8 a5 
1.5 FWHM (3.99 pixels), which was sufficiently large to = 

include most flux of a point source (more than 99.9% in our © 0.2 
case) and was not too large to include too many noisy empty T oo 


pixels. Most astrophysical studies require accurate flux 
measurement. It is important to check the reliability of the 
flux measurement. 

We used the fixed segment, i.e., the segment of the object 
(either galaxy or star) was fixed to the model image segment 
mentioned in previous subsections. If we use free segment or 
use resampled image, the conclusion does not change. We also 
note that Poisson resampling of the input flux was not used in 
our simulated images. In this way we make sure the difference 
only comes from the bias noise fluctuation within the same 
area. We will make a direct comparison between different noise 


Figure 8. Similar to Figure 5 but for row noise corrected images. 


types at the same nominal “readout noise” level and “signal-to- value. aie apathy ie = epy ihe aca 
foise ratio? will be normalized by the “true value” and presented in 
percentage. 
Figure 14 shows the measurement error as a function of 
3.6. Results 


average pixel flux in segment, which is defined as the total flux 
(in ADU) within the detection segment divided by the number 
of pixels of the segment. Here the position offset is defined as 
the absolute position offset with the sign of Y offset (i.e., 


We define the measurement error as the measurement 
difference between simulated image (model + noise) and pure 
model image, i.e., we treat the pure model result as “true 
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Figure 9. Similar to Figure 6 but for row noise corrected images. 


Table 1 
Random Parameters for Sérsic Models 
Parameter Min Max 
amplitude 100 10° 
effective radius (pixel) 3 10 
Sérsic index 1 5 
ellipticity 0 0.8 
position angle (degree) 0 360 
center x offset (pixel) —0.5 0.5 


center y offset (pixel) 0.5 L5 


|Ay|/AY x JAX? + AY?). Each data point is binned from 
1000 measurements. Error bars show the variation within the 
bin. Statistically, the mean value is accurate even for RN. But 
at any input brightness, the uncertainty of RN measurements 
are clearly worse than the Gaussian counterparts (note they 
have the same “readout noise”). Real galaxies and Sérsic 
models have similar results. 

Figure 15 shows the distribution of measurement error of all 
real galaxy images for six parameters. It shows the performance 
of RN is the worst for all parameters. The measurement error is 
especially large for Y position, segment flux and Kron flux. 


After correction, the performance of RN-corr is generally 
between Gauss-high and Gauss-low cases. It suggests that the 
row noise correction is not perfect but still very useful. 

The distribution is sometimes not symmetric. It may be 
related to the underlying algorithm and properties of the input 
galaxy templates. Such asymmetry does not affect our 
conclusion. Interestingly, the Y position uncertainty in RN 
case is significantly worse than X position. It is somehow 
expected, given the anisotropic nature of row noise. 

To better quantify the difference between different noise 
types, we plot the measurement uncertainty (the standard 
deviation of measurement error within the bin) as a function of 
average pixel flux in segment. Figures 16—18 show the results 
of all four noise types combined with three different models. 

The result is similar. Except the position angle of Sérsic 
images, RN is significantly worse than all other cases. RN-corr 
is much improved, comparable with Gaussian (between the 
high and low noise cases). The effect is much prominent in 
galaxies (including real galaxies and Sérsic models) than stars. 
It is probably due to the fact that the flux measurement 
apertures of stars are usually much smaller than that of 
galaxies. More pixels mean higher noise contribution from bias 
noise, and hence more vulnerable to row noise. It also clearly 
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Figure 10. The same object with different noise residual as background (displayed in the same color scale). Left column: real galaxy model; middle column: Sérsic 
model; right column: star model. The images are zoomed in for better visibility. 
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Figure 11. Real galaxy models used in simulation. Right ascension and decl. of the image center are shown for each galaxy (some galaxies are not centered). 
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Figure 12. Example Sérsic models used in simulation. Blue line shows the segment which contains 90% of total flux. 


demonstrates that the correction is useful but the corrected 
result is still not as good as pure white noise. 

By comparing the flux uncertainties with the same input 
brightness, it is found that the uncertainty of RN case is about 3-10 
times larger than Gauss-high for galaxies, and about 2 times larger 


for stars. In this paper, our test strictly controls the variables, leaving 
only the bias noise as the only parameter. So the effect of row noise 
we measured here can be directly translated into the increase of 
effective readout noise, by a factor of 2-10, depending on the target 
type. If we compare RN with Gauss-low, the difference is larger. 
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Figure 13. Example Gaussian models used in simulation. All panels are in the same scale. The images are zoomed in. 


4. Discussion and Summary 


In this paper, we studied the noise properties of bias images 
taken from an HR9090BSI sCMOS detector. We found a 
strong pixel cross-correlation in each row, i.e., the row noise, 
and some minor fluctuations with fixed vertical spatial 


frequencies in the column direction. A simple row noise 
correction was applied to the real bias image, leading to a much 
lower readout noise and much smaller pixel-to-pixel correla- 
tion. We made a set of simulated galaxies and stars with 
different noise residual images as background. A series of 
photometry tests were performed. The results suggest that the 
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Figure 14. Measurement error as a function of average pixel flux in detection segment (see text) for simulated stars. The blue data points are slightly shifted along X- 


axis for better visibility. 


row noise feature can significantly deteriorate the overall 
photometry performance, especially for galaxies. It effectively 
increases the readout noise by a factor of 2-10, depending on 
the type of object. 

We demonstrate that the existence of row noise pattern, i.e., 
pixel-to-pixel correlation, could greatly increase the photometric 
uncertainty. Such an effect is not simply due to the increase of 
nominal readout noise. In fact, at the same nominal readout noise 
level, the row noise case performs much worse than its Gaussian 
noise counterpart. It again proves the importance of pixel-to-pixel 
independence in astronomical high-precision photometry. It also 
makes a clear warning to anybody working on observation 
planning or camera design. The nominal “readout noise” may be 
misleading when the detector shows strong pixel-to-pixel 
correlation (no matter what the origin of such a correlation is). 
It may underestimate the required observation time for a given 
photometric accuracy, leading to an ambiguous observation 
result, or overestimate the performance of camera, leading to a 
failure of achieving the original design goal. 

We also note that, given the fact that the row noise is very 
common in sCMOS detectors, it is necessary to seriously 
consider the destriping method, either implemented in hard- 
ware or adopted at the data processing stage with efficient and 
robust algorithm. Our result suggests even the simplest 
destriping algorithm can make prominent improvement. Any 
future astronomical project using sCMOS detector should be 
aware of these facts. 

Our analysis also has limitations, although we note that some 
of them are out of the scope of this paper. First, our analysis of 
row noise is totally phenomenal. We treat the camera as a black 
box (due to some reasons). In order to fully understand the 
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origin of row noise in this specific CMOS model, we will need 
more knowledge of the detector system, including the 
manufacturing details of the chip and the readout electronics. 
Second, we only consider the additive signal-independent row 
noise. Other multiplicative and/or signal-dependent noise may 
also affect the photometry. All these additional noise 
components will further worsen the photometry result if not 
properly corrected. In order identify these noise components, 
we will need more data, especially images from real 
astronomical observations. We note that in near future the 
“Earth 2.0” project (Song et al. 2022) may give us a good 
chance. Third, this paper is focused on bias noise only. Our 
result is a good quantitative evaluator of the image readout 
performance. However, as a side effect, it makes the result 
looks too “exaggerated” when comparing to real situation. In 
practice, the bias noise usually only contributes a small fraction 
of the total photometric uncertainty. Many noise sources, such 
as Poisson noise, dark current, sky background and stray light, 
will make major contributions. Readers should be cautious 
when interpreting our results. To make a more realistic analysis 
of the photometric uncertainties, data from a more sophisticated 
simulation (e.g., Fang et al, in preparation, specifically 
designed for CSST'') should be used. Finally, we do not 
discuss too much about the algorithm for row noise correction. 
When more and more CMOS data appears in the field of optical 
astronomy, it may be interesting to systematically study the 
effect of existing methods to high-precision photometry based 
on real observation images. 


11 CSST simulation v2.0 release: https: //csst-tb.bao.ac.cn/code/csst_sim/ 
csst-simulation /-/tree/release_v2.0. 
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Figure 15. The distribution of measurement errors for real galaxy models. Different bias noise types are displayed in different colors. 
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Figure 16. The measurement uncertainty as a function of average pixel flux in detection segment for real galaxy models. Different bias noise types are displayed in 
different colors. 
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